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We present a theoretical estimate for a new observable: the cross-correlation between the Lyman-a 
flux fluctuations in quasar (QSO) spectra and the convergence of the cosmic microwave background 
(CMB) as measured along the same line-of-sight. As a first step toward the assessment of its 
detectability, we estimate the signal-to-noise ratio using linear theory. Although the signal-to-noise 
is small for a single line-of-sight and peaks at somewhat smaller redshifts than those probed by the 
Lyman-a forest, we estimate a total signal-to-noise of 9 for cross-correlating QSO spectra of SDSS- 
III with Planck and 20 for cross-correlating with a future polarization based CMB experiment. The 
detection of this effect would be a direct measure of the neutral hydrogen-matter cross-correlation 
and could provide important information on the growth of structures at large scales in a redshift 
range which is still poorly probed by observations. 
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Introduction. The Lyman-a forest - the absorption 
seen in quasar (QSO) spectra caused by intervening neu- 
tral hydrogen of the intergalactic medium (IGM) - has 
the potential to provide precise information about the 
matter distribution down to small scales. In the near 
future, the SDSS III Baryon Oscillation Spectroscopic 
Survey (BOSS) will measure absorption spectra toward 
160, 000 QSOs, a fiftyfold improvement on existing sur- 
veys This large sample could potentially provide 
unprecedented constraints on neutrino masses 
plitudc and slope of the matter power spectrum 
flationary parameters and the running spectral index 
Furthermore, the baryon oscillations in the Lyman-a for- 
est can be used to constrain dark energy and curvature 
and high-resolution QSO spectra can also play a sig- 
nificant role in constraining the properties of dark matter 
on small scales @. 

The theoretical studies assume that the physics of the 
Lyman-a forest is relatively simple and that the observed 
flux fluctuations faithfully trace the dark matter distribu- 
tion. Hydrodynamical simulations play an essential role 
in testing this assumption through comparisons with data 
and in calibrating the measurements. However, many 
systematic errors still pose challenges at the sub-percent 
level of accuracy required by the data. The main sources 
of uncertainty include, but are not limited to: the ac- 
curate modeling of fluctuations in the ionization back- 
ground, the QSO continuum fluctuations, the uncertain- 
ties in the slope of the temperature-density relation in 
the IGM and fluctuations about this mean relation, the 



effect of galactic supcrwinds on the Lyman-a flux and 
the contamination of the Lyman-a spectrum with metal 
lines (e.g. Q). The ingredients that go into the simula- 
tions could be better controlled if we could directly and 
independently test how the fluctuations in the transmit- 
ted flux and the dark matter are related. 

In this Letter, we propose one such observable: the 
cross-correlation of the Lyman-a flux fluctuations with 
the convergence field extracted from the gravitational 
lensing of the cosmic microwave background (CMB). 
Gravitational lensing of the CMB is caused by the deflec- 
tion of CMB photons by intervening large scale structure 
Q. Lensing breaks the statistical isotropy of the CMB, 
introduces non-Gaussianitics and produces B-mode po- 
larization from E- modes. Observations of small scale 
temperature and polarization fluctuations can be used to 
reconstruct the effective deflection field that lensed the 
CMB WMAP measurements have enabled the first 
detection of CMB lensing [l(J 11 1. With the ongoing 



and upcoming high resolution temperature and polariza- 
tion based CMB experiments like the Atacama Cosmol- 
ogy Telescope, the South Pole Telescope, Planck, QUIET 
and PolarBearR [13], high fidelity lensing reconstruction 
will soon become a reality. The timely convergence of the 
Lyman-a surveys and high resolution CMB experiments 
in coming few years will make the measurement of their 
cross-correlation a possibility. 

We note that cross-correlating CMB temperature 
maps with the large-scale structure (LSS) is currently 
done with many LSS probes including Lyman-a (e.g. 
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[l3(), however we focus here on the convergence of the 
CMB rather than the CMB temperature fluctuations. 
The cross-correlation signal between these two quanti- 
ties should provide insights on the relative bias between 
matter and flux (in a way somewhat analogous to [Hj ) 
and may become a powerful tool to calibrate the rela- 
tionship between the IGM and the dark matter. As a 
first step toward understanding the information content 
and the detectability of this correlation, we provide here 
a simple theoretical estimate of the signal and of the sig- 
nificance with which it can be measured by the upcoming 
experiments. 

Evaluation of the correlators. We calculate the corre- 
lations between the fluctuations in the flux 8J-(n, x) an d 
their variance 8T 2 (h, \) measured for a single QSO - 
whose line-of-sight (los) lies in the direction n and whose 
spectrum covers a redshift range Az - and the CMB con- 
vergence n(n) measured in the same direction. While the 
convergence measures the projected matter fluctuations 
along the los, the quantity 8T is sensitive to the typi- 
cal fluctuations of the Lyman-a forest which are in the 
range 1-100 comoving ft, _1 Mpc. Since structures grow 
at a faster rate in overdense regions (which in turn are 
characterized by a positive convergence) the cross- 
correlation of the convergence with the flux variance is 
expected to increase the signal. Also, using the flux vari- 
ance makes the observable less sensitive to the continuum 
fitting uncertainties of the QSO spectrum since 8T re- 
quires correct extrapolation of the continuum from long- 
ward of Lyman-a . 

The effective CMB convergence measured in the di- 
rection n depends on the amplitude of matter density 
fluctuations S(n,x) along the los and reads: 
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where the integral is along the los, xls is the comov- 
ing distance to the last scattering surface, Wl{ X ,Xls) = 
x{xls ~x)/xls is the usual lensing window function and 
a(x) is the scale factor at comoving distance x- 

We assume that on large scales fluctuations in the opti- 
cal depth are linearly proportional to fluctuations in the 
matter density [l5j], motivated by the so-called fluctu- 
ating Gunn-Peterson approximation [l6l . 17 1 that relates 
flux and matter F = exp [— A(l + S) 13 ], where A and f3 are 
two redshift dependent functions: A is of order unity and 
is related to the mean flux level, baryon fraction, IGM 
temperature, cosmological parameters and the photoion- 
ization rate of hydrogen, while f3 depends on to the so- 
called IGM temperature-density relation (e.g. 0, UK]). 
However, for the purpose of this Letter, it is fair to neglect 
their evolution with redshift and to assume that they are 
constant over the observed QSO redshift interval in each 
los. In the large scale limit, the exponent is small and we 
can expand it. Keeping only up to the lowest order term 
in 5 [l7|, we get: 8T(n,x) » -A(38(h,x)- We stress 



that this approximation is valid in linear theory neglect- 
ing not only the non-linearities produced by gravitational 
instabilities but also those introduced by the definition 
of the flux, thermal broadening and peculiar velocities. 
Higher order terms in the expansion will, in general, con- 
tribute to the cross-correlations being estimated here. To 
keep the calculations analytically tractable, in this Letter 
we neglect these higher order terms. 

We define the m— th moment of the flux fluctuation 
averaged over the los interval probed by the Lyman-a 
spectrum (extending from X i to xq) as: 
rXQ 

8T m {n) = / d X 6T m (n, X ) 

■>Xi 

rXQ 

« (-Af3) m / d x 6 m (n, X ). (2) 

The correlation of this quantity with the CMB conver- 
gence k measured along the same los is then, 
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We focus here on the case m = 1 (and m = 2), which 
is the cross-correlation between the CMB convergence 
and the QSO mean flux fluctuation (and variance). To 
proceed further we Fourier transform the cumulant corre- 
lator (8 m (n, x q ) 8(n, Xc)}- The case m = 1 reduces to the 
2-point correlation function of the dark matter overdensi- 
ties along the los, while the case m = 2 requires the use of 
the 3-point correlation function or of its Fourier counter- 
part, the bispectrum [20l |. In both cases we also include 
in the analysis the fact that both the Lyman-a spectra 
and the CMB convergence arc being measured with a 
finite resolution. To this end we add appropriate gaus- 
sian window functions that limit the modes contributing 
to the CMB convergence (perpendicular to the los) to 
wavenumbers \k±\ < kc and the ones contributing to the 
Lyman-a spectra (parallel to the los) to fc; < fcii < kj,. 
When these cutoffs are included, the treatment of the 
m — 2 case becomes involved. However, by expanding in 
power series the modified Bessel functions arising from 
the angular integration in fc-space, it is possible to ob- 
tain an exact series solution for the cumulant correlator. 

We calculate the value of (8J- k) and of (8 J- 2 n) as 
a function of the QSO redshift z and of the redshift 
range Az spanned by the Lyman-a spectrum for a flat 
ACDM cosmology with fl 0m = 0.25, <j 8 = 0.84, n s = 0.96 
and Hq = 72 km/s/Mpc. The results obtained - shown 
in Fig. [I]- make physical sense: if the Lyman-a spectrum 
covers a larger redshift range the correlators increase, as 
more information is carried by the longer spectrum. On 
the other hand, if the source QSO redshift is increased 
while keeping the length of the spectrum fixed the cor- 
relators become smaller because the Lyman-a spectrum 
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FIG. 1: Values of the correlators {8J-k) (left panel) and (8T 2 k) (right panel) along a single los as a function of the source redshift 
z and of the length of the measured spectrum Az. The cutoff scales assumed are ki = 0.01 hMpc- 1 and k L = 4.8/iMpc" 1 for 
the QSO spectrum and kc = 0.021 /iMpc -1 for the CMB convergence, consistent with SDSS-III and Planck specifications. We 
set A = /3 = 1 so that it is straightforward to rescale the present figures to match the measured values of such parameters. 



is probing regions where structure had less time to form. 
Furthermore, the rate at which the correlators are in- 
creasing as a function of Az is itself increasing with 
the CMB and Lyman-a experiment resolutions fej, and 
kc- Finally, as expected the (SJ^k) correlator is about 
two orders of magnitude larger than (SJ-k): regions of 
higher convergence have typically higher mean density 
and hence larger fluctuations' variance due to enhanced 
growth [13]. 

Estimate for the correlators' signal-to-noise ratio. To 
assess the detectability of the above correlations by 
CMB experiments and QSO surveys we need to evalu- 
ate their signal-to-noise ratio (S/N). The calculation of 
the variances of ST K and ST 2 k requires the evaluation 
of (5 J- 2 k 2 ) and of (8!F a k 2 ), which are proportional to 
(SqSgi S C 5 C >) and {5 2 5 2 ,5 C S C >), respectively [we adopt the 
notation Si = 5(xi)]. The exact evaluation of the lat- 
ter requires the calculation of a six-point function, the 
form of which to our knowledge has never been obtained. 
However, an estimate for the dominant contributions to 
the variance can be obtained using the approximation 

(5 2 5 2 ,5 C S C ,) « 2{5 2 S C ) 2 + {S C S C ,) ({5 2 ){S 2 q ,) + 2{S q 5 q ,) 2 ) ,(4) 

The estimates for the S/N ratio for the measurement 
of ST k and 5 J- 2 k along a single los are shown in Fig. [2] 
As before, increasing the length of the spectrum or de- 
creasing the redshift of the source will increase the S/N 
ratio, as more information is added by the longer spec- 
trum in the first case and as the shallower redshift will 
allow to probe regions with clumpier structure in the sec- 
ond case. Also, comparing the two panels it is possible to 
notice how the correlators' different dependence on the 
growth of structure affects the redshift range over which 
a given S/N ratio can be obtained. Finally, we also no- 
tice that the S/N for ST k is about four times larger than 



the one for b~T 2 k, a consequence of the fact that cosmic 
variance is generally higher for higher order statistics. 

Conclusions. We presented a calculation of the cross- 
correlation signal (and its noise) between the (variance 
of the) Lyman-a flux fluctuations and the CMB conver- 
gence. We found that the cross-correlation signal peaks 
at somewhat smaller redshifts than those usually probed 
by Lyman-a alone, but it will be detectable with future 
data sets with a large number of QSO. In particular, we 
estimate that the S/N ratio for (8 J- 2 k) obtained by cross- 
correlating the BOSS QSO sample with Planck data will 
yield a S/N ratio of 9.6. The use of higher definition data 
arising from a hypothetical CMB polarization experiment 
covering 8000 sq. degrees would further increase the es- 
timated S/N to 20. In all the above cases the present 
analysis shows that this correlation will be observable. 

A few caveats are also in order. First, as large con- 
tributions to the signal may actually be arising at small 
separations, it is necessary to stress that the estimates 
obtained for the variance of the correlators - and con- 
sequently for S/N - do not take into account the non- 
linearities characteristic of small separations. The latter 
however are not expected to alter significantly the results, 
given the QSO mean redshift. Second, it is necessary to 
point out that the estimates of Fig.[5]do not take into ac- 
count the detector systcmatics, which in a more realistic 
calculation should also be included. Finally, we note that 
the estimators considered may not be the ones character- 
ized by the least variance. It is reasonable to speculate 
that the use of redshift weighted estimators in conjuction 
with lower redshift convergence maps obtained through 
LSS surveys may allow an improvement in the S/N ra- 
tio. The formalism developed can further be extended to 
cross-correlate the Lyman-a flux with convergence maps 
obtained from optical lensing surveys once these will be- 
come available at the redshifts of the Lyman-a forest. 
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FIG. 2: Estimate of the S/N ratio for the measurement of the correlators (STn) (left panel) and {8T 2 k) (right panel) along a 
single los as a function of the source redshift z and of the length of the measured spectrum Az. The cutoff scales assumed are 
ki = 0.01 ftMpc -1 and kt = 4.8/iMpc -1 for the QSO spectrum and kc = 0.021 /iMpc -1 for the CMB convergence, consistent 
with SDSS-III and Planck specifications. 



This would result in an improvement in S/N ratio. 

Finally, the present work is just a first step in assessing 
the information content and the detectability of Lyman- 
a -CMB lensing correlations and more refined calcula- 
tions, based on non-linear theory and hydrodynamical 
simulations should be performed to address this more 
quantitatively, as future observational programs should 
enable to perform such a measurement with several con- 
sequences. First, this would provide independent con- 
straints on the bias between flux and matter, additional 
to the ones that are obtained by high resolution QSO 
spectra and hydrodynamical simulations. Second, being 
sensitive to the growth of structures at intermediate to 
small scales, this correlation could be used to test for 
models of early dark energy plj or modification of grav- 
ity that produce scale dependent growth evolutions Q 
and to provide constraints on the neutrino masses. 
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